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Abstract 

The  dynamic  performance  of  PEM  fuel  cell  is  one  of  the  most  important  criteria  in  its  design  with  application  to  mobile  systems  and  portable 
devices.  To  analyze  the  features,  this  work  employs  an  unsteady  model  about  single  phase  transport  in  cathode  side  of  PEM  fuel  cell  with 
interdigitated  flow  field,  which  considers  both  convection  and  diffusion  processes.  Two  types  of  dynamic  performances,  start-up  and  state- to- state 
operations,  are  analyzed.  The  effects  of  channel  width  fraction,  porosity  of  the  gas  diffusion  layer,  pressure  drop  and  the  surface  overpotential  of 
the  catalyst  layer  on  the  dynamic  performance  are  investigated  in  detail.  Predicted  results  found  that  the  response  time  is  generally  quite  fast,  less 
than  0.1  s,  to  reach  the  90%  response.  The  time  interval  from  the  start-up  to  a  steady  state  or  from  one  steady  state  to  another  steady  state  mainly 
depends  on  the  end  condition. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Proton  exchange  membrane  (PEM)  fuel  cells  have  extensive 
applications  due  to  their  merits  of  high  efficiency,  low  pollu¬ 
tant  emission  and  convenience,  and  now  the  main  application 
is  as  power  generators  of  mobile  systems.  Moreover,  it  also 
has  some  promising  applications  in  portable  devices,  such  as 
mobile  phones,  laptop  computers  and  digital  cameras.  In  appli¬ 
cations  of  PEM  fuel  cells  to  these  systems,  the  dynamic  response 
of  the  system  is  important.  The  dynamic  response  involves 
two  operations,  start-up  operation  and  state-to- state  operation. 
The  time  interval  from  the  start  to  a  steady  state  or  from  one 
steady  state  to  another  steady  state  is  one  of  the  significant  per¬ 
formance  characteristics  of  the  system.  The  dynamic  response 
of  the  fuel  cell  mainly  depends  on  operating  conditions,  mass 
transport  process  in  the  cells,  electrochemical  reaction  kinetics, 
mechanical  design  and  manufacturing  process,  etc.  Among  these 
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factors,  the  reactant  gas  transport  in  the  cell  is  one  of  the  crucial 
points. 

The  gas  diffusion  layer  (GDL)  morphology,  geometry  and 
porosity,  as  well  as  the  flow  channel  arrangement,  e.g.,  parallel, 
serpentine  or  interdigitated  channel,  gas  pressure  drop  have  sig¬ 
nificant  influences  on  the  optimization  of  the  cell  performance. 
Some  studies  on  modeling  mass  transport  and  parametric  anal¬ 
ysis  in  the  PEM  fuel  cell  have  been  performed.  The  pioneer 
work  is  from  Bernardi  and  Verbrugge  [1],  who  developed  a  one¬ 
dimensional  isothermal  model  to  examine  mass  transport  in  fuel 
cell.  Gurau  et  al.  [2]  proposed  a  one-dimensional  model  to  exam¬ 
ine  the  characteristics  of  the  mass  transport  of  reactant  gas  in  a 
half-cell,  in  which  the  effects  of  the  porosity  and  the  tortuosity  of 
GDL  and  the  catalyst  layer  were  investigated.  Springer  et  al.  [3] 
predicted  the  net  water-per-proton  flux  ratio  across  the  mem¬ 
brane.  Okada  and  colleagues  [4]  studied  water  transport  with 
and  without  impurity  ions  in  membranes  of  PEM  fuel  cells.  For 
more  details  of  the  steady-state  flow  distribution  and  gas  diffu¬ 
sion,  a  number  of  studies  on  two  and  three  dimensional  transport 
phenomena  in  PEM  fuel  cells  have  been  carried  out.  Dutta  et  al. 
[5]  used  the  software  package  to  analyze  the  distributions  of  gas 


J.  Zou  et  al.  /  Journal  of  Power  Sources  159  (2006)  514-523 


515 


Nomenclature 

D 

mass  diffusivity,  m2  s-1 

DQ 

effective  mass  diffusivity,  m2  s-1 

F 

Faraday  constant,  96487  C  mol-1 

j 

transfer  current  density,  Am-2 

i 

average  current  density,  Am-2 

k 

permeability,  m2 

k 

length  of  channel,  m 

k 

length  of  shoulder,  m 

Mo2 

molecular  mass  of  oxygen,  32.0kgkmol_1 

Mu2o 

molecular  mass  of  oxygen,  18.0kgkmol_1 

n 

number  of  electrons  in  relation  participating 

P 

pressure,  Pa 

R 

universal  gas  constant,  8.314  J  mol-1  K-1 

hydaulic  radius,  m 

Sc 

thickness  of  catalyst  layer,  m 

Sd 

thickness  of  gas  diffusion  layer,  m 

sm 

thickness  of  membrane,  m 

t 

time,  s 

tc 

time  constant 

T 

temperature,  K 

Vx 

superficial  velocity  in  v  direction 

Vy 

superficial  velocity  in  y  direction 

V 

cell  voltage 

Voc 

ideal  cell  potential,  V 

Wo2 

mass  fraction  of  oxygen 

Wu20 

mass  fraction  of  water  vapor 

X 

v  direction 

y 

y  direction 

Greek  letters 

o\ 

porosity  of  GDL 

02 

porosity  of  catalyst  layer 

T] 

surface  overpotential,  V 

A 

channel  width  fraction,  lc/ls 

Aw 

membrane  water  content 

P 

Viscosity 

K 

reaction  rate  constant,  molm-3  s-1 

P 

density,  Kg  m-3 

O 

ionic  conductivity  of  the  ionomer,  £2_1  m-1 

X 

tortuosity 

0 

ionomer  phase  potential,  V 

density  and  velocity  in  flow  channel  and  GDL.  Jordan  et  al.  [6] 
examined  the  influence  of  diffusion-layer  morphology  on  cell 
performance.  A  model  of  the  hydrophobilicity  and  the  porosity 
of  the  gas-diffusion  layer  was  used  to  explain  the  influence  of 
the  diffusion-layer  morphology.  Yi  and  Nguyen  [7]  developed 
a  two-dimension  model  to  examine  the  multi-component  gas 
transport  in  PEM  fuel  cells  of  interdigitated  flow  channel.  And 
then  they  [8]  developed  a  two-phase  model  to  examine  the  water 
and  thermal  management  of  PEM  fuel  cells  with  interdigitated 
flow  field.  They  found  that  the  interdigitated  gas  distributors 
can  reduce  mass-transport  overpotential  and  improve  the  flood¬ 


ing  phenomena  in  the  cathode  side.  Yan  et  al.  [9]  developed  a 
two-dimensional  model  to  investigate  gas  reactant  transport  at 
various  conditions  of  flow  channel  width  ratio  (A  =  channel  width 
/c/shoulder  width  /s)  and  GDL  porosity  £\ .  They  disclosed  that 
an  increase  in  either  A  or  £i  may  lead  to  a  better  cell  perfor¬ 
mance.  At  the  case  with  a  relatively  low  overpotential,  better 
uniformity  in  current  density  distribution  along  the  width  of 
the  cell  can  be  attained.  Chen  et  al.  [10]  investigated  the  tran¬ 
sient  behavior  of  the  water  transport  across  the  membrane  of  the 
PEM  fuel  cell  to  seek  for  effective  control  schemes  so  that  the 
best  dynamic  performance  of  the  fuel  cell  can  be  obtained.  It  is 
found  that  both  a  larger  starting  operational  current  density  and  a 
smaller  operational  current  density  can  lead  to  a  shorter  dynamic 
response  time,  while  control  of  the  humidification  of  the  fed  fuel 
is  the  most  powerful  as  well  as  the  most  feasible  scheme.  Yan 
et  al.  [11]  studied  transient  gas  transport  in  the  cathode  side 
of  PEM  fuel  cell  with  straight  flow  field.  They  concluded  that 
an  increase  in  channel  width  ratio  A,  or  surface  overpotential  ij 
may  lead  to  a  faster  dynamic  response  of  a  PEM  fuel  cell  sub¬ 
jected  a  start-up  operation,  and  a  higher  surface  overpotential 
results  in  a  fuel  cell  with  a  faster  surface  reaction  rate.  The  time 
interval  from  start-up  to  a  steady  state  may  be  as  long  as  10  s 
due  to  the  mass  transport  lag,  but  it  is  only  about  0.4  s  to  reach 
the  90%  response.  The  surface  overpotential  ij  represents  the 
voltage  loss  due  to  chemical  reaction  within  the  catalyst  layer, 
and  is  the  effect  of  the  current  of  a  fuel  cell.  A  higher  surface 
overpotential  means  a  faster  surface  reaction  rate  when  other 
conditions  are  the  same.  For  a  certain  fuel  cell  system,  the  value 
of  r)  or  the  reaction  rate  can  be  controlled  by  adjusting  external 
load. 

Although  the  previous  studies  concern  many  issues  about  the 
PEM  fuel  cells,  from  steady  state  analysis  to  unsteady  analysis, 
no  work  has  been  done  about  dynamic  analysis  of  fuel  cell  with 
interdigitated  flow  field.  It  features  the  combined  influences  of 
convection  and  diffusion  processes  and  excellent  anti-flooding 
function.  The  objective  of  this  work  is  to  numerically  inves¬ 
tigate  the  dynamical  response  of  this  type  of  PEM  fuel  cells. 
Parametric  analysis  is  performed.  The  emphasis  of  the  analy¬ 
sis  is  placed  on  the  influences  of  the  channel  width  ratio  A,  the 
porosity  of  the  GDL  s\  and  pressure  drop  A P  on  fuel  cell.  By 
applying  the  time-dependent  mass  transport  model  developed 
in  this  work,  the  effects  of  these  parameters  and  the  surface 
overpotential  r\  on  the  transient  characteristics  of  the  mass  trans¬ 
port  and  the  dynamic  response  of  the  cell  performance  are 
explored. 

2.  Numerical  modeling 

2.7.  Problem  statement 

In  this  work,  the  single  phase  multi-component  transport  pro¬ 
cess  is  considered  only.  This  assumption  has  its  limitation.  The 
maximum  vapor  pressure  in  the  cell  is  limited  to  be  below  the  sat¬ 
urated  pressure  at  the  operating  temperature,  which  is  0.047  MPa 
at  80  °C.  And  generally,  the  value  of  the  maximum  vapor  pres¬ 
sure  increases  with  current  density  when  it  is  below  the  saturated 
pressure.  Therefore,  the  density  current  can  not  be  too  large. 
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Fig.  1.  Half-cell  model  of  a  PEM  fuel  cell.  BP:  bipolar  plate;  GDL:  gas  diffusion  layer;  CL:  catalyst  layer;  /c:  width  of  the  channel;  /s:  width  of  the  shoulder. 


Otherwise,  the  result  will  not  reflect  the  real  flow  process.  Wang 
et  al.  [12]  predicted  the  onset  current  density  for  the  two-phase 
regime  at  different  inlet  dry  air  velocity  under  certain  configu¬ 
ration.  When  the  velocity  is  about  1ms-1,  the  critical  current 
density  is  about  2  A  cm-2.  Although  single  phase  condition  is 
not  correct  at  the  operating  condition  of  high  current  density, 
it  is  often  reasonable  in  low  current  density  with  dry  air  entry 
when  no  liquid  water  is  formed  or  the  liquid  saturation  is  very 
low. 

Multicomponent  consideration  accounts  for  the  effects  of 
such  inert  species  as  nitrogen  and  water  vapor.  Two-dimensional 
condition  is  assumed,  and  the  computational  domain  is  perpen¬ 
dicular  to  the  axial  direction  of  channels  in  bipolar  plate.  This 
treatment  is  consistent  with  the  main  mass  transport  direction 
normal  to  the  channels  in  the  fuel  cell.  As  shown  in  Fig.  1,  only 
half-cell  is  adopted.  Gas  mixture  enters  the  gas  diffusion  layer 
(GDL)  from  channels  connected  with  inlet  plenum  of  bipolar 
plate  and  exits  from  its  neighbor  channels.  For  simplicity,  a  unit 
structure  is  considered  as  the  computational  domain,  as  shown 
in  the  right  graph.  There  are  totally  three  layers,  gas  diffusion 
layer  (GDL),  catalyst  layer  (CL)  and  membrane.  Gas  transports 
within  the  GDL  and  CL,  and  convection  and  diffusion  processes 
are  both  crucial  in  this  problem,  which  is  different  from  the 
gas  transport  in  the  PEM  fuel  cell  with  other  types  of  flow 
fields. 

2.2.  Assumption 


2.3.1.  Continuity  equation 
d(PVx)  _  d(pvy)  _  j 

dx  3  y  2  nF 


(Mq22Mh2o)  +  £ 


dp 

~dtm 


Where  s  represents  the  porosity  of  GDL  £i  or  CL  s 2 ,  and  j 
represents  the  cathode  transfer  volumetric  current  density  and 
can  be  formulated  by  the  following  equation  according  to  the 
ORR  assumption. 


nFKWo2  exp 


0 


in  the  catalyst  layer 
in  the  GDL 


2.3.2.  Momentum  equation 

According  to  Darcy’s  law,  superficial  velocity  in  porous 
media  is  expressed  as  below. 


k  3  P 
pi  dx 

k  3  P 
pi  3  y  ' 


(3a) 

(3b) 


In  Eqs.  (3),  the  k  is  permeability  of  porous  media.  To  test 
the  influence  of  porosity  of  GDL  on  mass  transfer,  especially 
convection  process,  Carman-Kozeny  equation  [13]  is  employed 
to  express  the  permeability  k  as  a  function  of  porosity  s  and 
hydraulic  radius  R/t. 


In  this  study,  the  membrane,  GDL,  and  CL  are  respectively 
treated  as  homogeneous  porous  media  with  uniform  morpho¬ 
logical  properties,  such  as  porosity,  tortuosity  and  permeability. 
The  gas  mixture  is  considered  to  be  ideal  gas  and  the  viscosity  is 
assumed  constant.  Isothermal  condition  is  assumed,  so  energy 
equation  can  be  omitted. 

Gravity  effect  is  neglected.  Darcy’s  law  is  employed  to 
account  for  superficial  velocity  in  porous  medium.  Fick- 
ian  law  is  introduced  to  account  for  species  diffusion  pro¬ 
cess.  The  oxygen  reduction  reaction  (ORR)  is  adopted  in  the 
analysis. 

2.3.  Mathematical  model 


where  the /is  the  Kozeny  constant.  It  is  assumed  that  the  mor¬ 
phology  of  the  same  kind  of  porous  media  keeps  similar  when 
its  porosity  changes,  and 


Rh  =  C\fe 


So  the  relation  between  permeability  and  porosity  is  established 
as  below 


where  f 


f 

C 


The  single  phase  multicomponent  model  of  PEMFC  con¬ 
sists  of  four  kinds  of  principle  equations:  continuity,  momentum, 
species  and  phase  potential  equations. 


In  this  work,  the  value  of  /  is  set  to  be  3  x  10nm  2. 
When  the  value  of  s  is  0.6,  the  corresponding  permeability  is 
1.2  x  10-12  m2. 
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2.3.3.  Species  equation 

Only  two  kinds  of  species  equations  are  needed,  and  here 
oxygen  and  water  vapor  are  selected. 


9 

dx 


pD0 


e  ^w02 


dx 


+ 


9 

dy 


PDq 


9  j 

~  —(pwo2vy)  =  —  —  Mq2  +  e 
dy  2nF 


dwo2  \  9 

■— )  -  3i(PW0‘V’) 
d(pwo2) 


3 1 


(7) 


3 

dx 


Pdh2o 


dwn2o 

dx 


9 

+  jr(pD 


h2o 


9 


—  (PWH2OVy) 
dy 


=  £ 


dy 

d(pwu2o ) 

dt 


9u;h2o 

dy  7  dx 


9 

)  -  —(pwH2OVx) 


— 4mH2° 

nF 


(8) 


Combining  Eqs.  (1),  (3),  (7)  and  (8)  yields 


9 

dx 


pDQ 


e  dwo2 


+ 


9 


2  dx  J  dy 
dP 


pDQ 


e  dwo2 


dy 


+ 


k  dwo , 


pi  dx 


P 


dP 

dx 


k  dwoo 

H - I  p 

It  dy  \  dy 

+  2u;h2o^h2o)  +  ps 


J 


2  nF 
dwo2 


(Mo2(1  wo2 ) 


dt 


(9) 


k  dwu2o  f  dP\  k  dwu2o  f  dP\ 
p  dx  \  dx  J  /j  dy  \  9y  J 

j  ,  ,, 


+  /9£ 


2nF 
dwu2o 


(-Mq2wu2o  +  2u;h2o^h2o  2Mh2o) 


9t 


(10) 


The  effective  diffusivities  of  oxygen  and  water  vapor  are  as 
below. 


Dq2  —  Do2st 

dh2o  =  Dn2o£T 


2.3.4.  Phase  potential  equations 
9  /  d0\  9/90 

—  o  —  H - o  —  I  =  0  in  the  membrane 

dx  \  dx  J  dy  \  dy 


(Ha) 

(Hb) 


9  /  r  90 
so 


dx 


dx 


+ 


9 

dy 


os 


90 

dy 


=  j  in  the  catalyst  layer, 


(13) 


In  the  above  equations,  the  0  is  ionomer  phase  potential  and 
o  is  the  membrane  ionic  conductivity  and  is  evaluated  by 


o  =  (0.5139AW  —  0.326)  x  exp 


1268  x 


1 


303 


1 

T 


(14) 


2.3.5.  Boundary  and  initial  conditions 

For  description  simplicity,  each  boundary  is  numbered  in 
Fig.  1,  from  B1  to  B 8.  At  boundary  Bl, 


P  =  P{n,  U)q2  =  Wo2in,  W\\2o  =  ^H20,in 


(15a) 


In  this  work,  wo2m,  ^H2oin  are  respectively  set  0. 17  and  0.01 
At  boundary  B2, 


P=  R 


out- 


dwo2 

dn 


=  0, 


dwn 


20 


dn 


=  0 


(15b) 


Here  Pout  is  set  0. 1  MPa. 

At  boundary  B7,  total  mass  flux  is 


./  J 


k  dP 

—p —  =  2Sca 
p  dy  nF 


Mu2o 


(15c) 


where  ar  represents  the  net  water  transport  coefficient  in  the 
membrane,  and  is  mainly  combined  effect  of  electro-osmotic 
drag  and  diffusion.  The  value  of  a'  depends  on  several  condi¬ 
tions,  such  as  type  and  thickness  of  membrane,  water  content  of 
membrane,  and  have  a  quite  wide  range  [14, 15] .  For  certain  con¬ 
ditions,  the  value  of  a'  can  be  assumed  constant,  and  is  specified 
to  be  0.5  in  this  work  according  to  related  stusdies  [7,16]. 

Oxygen  and  water  vapor  flux  are  respectively 


pDQ 


e  dwo2 


dy 


k  dP 
+  w0l-p— 
p  dy 


=  0 


(15d) 


B7 


9u;h20  ,  k  dP 

/^H20  ^7  1“  wR20~P 


dy 


r  dy 


B7 


=  2Sca'  2-Mu2o 
nF 


(15e) 


where  Sc  is  height  of  the  catalyst  layer.  Here  there  is  an  assump¬ 
tion  that  electrical  conduction  is  one-dimensional  and  oxygen 
concentration  at  the  interface  represents  average  value  in  the 
catalyst  layer. 

Combining  Eqs.  (15c)-(15e),  we  have  the  reduced  expression 
for  the  condition  of  boundary  B7 


dP 

dn 


B7 


dwo: 

dn 


B7 


(12)  9u;h2o 
dn 


2Sca'(j/nF)Mu2o 
( k/m)p 

_  -2wo2Scaf(j/nF)Mu2o 

(1  -  u)ft2o)2ScaXj / nF)M\i2o 


(15f) 


(15g) 


B7 


pDQ 


(15h) 


h2o 


At  boundary  B6, 


pDQ 


e  dwo2 


dn 


=  pd  b 


e  dwo2 


+ 


2  dn 


dn 


+ 


k  3  P 


3  n 


k  3  P 
I  v/x  dn 


(15i) 


(15j) 


(15k) 
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The  boundary  conditions  for  phase  potential  0  are  as  follows. 
At  boundaries  B4  and  B5 


90 


=  0 


dx 

At  boundary  B7 


0+  =  0-, 


90 


dy 


90 


dy 

At  boundary  B6 
90 

=  0 


+ 


dy 

At  boundary  B8 
0  =  0 


(16a) 


(16b) 


(16c) 


(16d) 


(16e) 


There  are  two  types  of  initial  condition.  For  start-up  opera¬ 
tion,  the  initial  condition  is  as  below. 


wo2  —  0.0,  w  h2o  =  0-0  (17) 

And  for  state-to- state  operation,  the  initial  condition  is  the 
gas  distribution  of  the  former  steady  state. 

The  cell  potential  V  is  evaluated  by  subtracting  surface  over¬ 
potential  r)  and  the  average  phase  potential  difference  A0 
between  boundary  B  8  and  boundary  B7  from  the  ideal  cell  poten¬ 
tial  Voc,  that  is, 

V  =  Vqc  ~  0  ~  A0.  (18) 


3.  Numerical  method 

To  numerically  solve  the  above  set  of  equations,  the  implicit 
finite  difference  method  is  adopted.  The  time  derivative 
term  is  approximated  by  the  implicit  or  backward  difference 


scheme,  and  space  derivative  term  by  the  central  difference 
scheme.  With  the  coupled  finite-difference  equations,  the  linear 
algebraic  equations  are  solved  simultaneously,  firstly  obtaining 
approximate  pressure  and  velocity  profile  by  solving  continuity 
and  momentum  equations  and  then  calculating  species  profile 
by  species  equations.  The  uniform  grid  system  is  used  in  this 
study.  The  grid  points  in  the  x  and  y  directions  are  81  x  316, 
respectively,  which  is  tested  to  satisfy  the  grid  independent 
requirement. 

4.  Results  and  discussion 

Four  important  parameters  are  chosen  to  investigate  their 
effects  on  the  performance  of  the  PEM  fuel  cells.  They  are 
channel  width  fraction  A  =  /c//s,  the  porosity  of  GDL  s  \ ,  pres¬ 
sure  drop  A P  and  surface  overpotential  77.  The  values  of  other 
parameters  considered  in  the  present  study,  listed  in  Table  1 ,  are 
derived  from  the  practical  applications.  For  the  PEM  fuel  cell 
with  straight  flow  field,  channel  width  fraction  A  affects  the  cell 
performance  significantly,  and  an  increase  in  A  may  lead  to  a 
better  cell  performance  [9,1 1].  In  this  work,  however,  the  results 
show  that  X  just  slightly  influences  the  cell  performance.  The 
different  transport  process  and  boundary  conditions  may  lead  to 
this  result.  For  the  PEM  fuel  cell  with  straight  flow  field,  only 
diffusion  process  accounts  for  gas  transport  and  oxygen  mass 
fraction  wo2  at  the  exit  is  considered  to  be  constant.  While  for 
the  PEM  fuel  cell  with  interdigitated  flow  field,  the  convection 
and  the  diffusion  processes  are  combined  and  the  normal  gradi¬ 
ent  of  wo2  at  the  exit  is  zero.  So,  this  paper  will  not  analyze  the 
effect  of  A  further  on  the  cell  performance,  and  the  default  value 
of  A  is  fixed.  The  results  are  organized  in  two  sections,  start-up 
operation  and  state-to- state  operation. 

Just  as  Yan  et  al.  [11]  indicated,  although  it  may  spend  con¬ 
siderably  long  time  to  reach  completely  a  steady  state  due  to 
the  mass  transport  lag,  it  needs  only  quite  little  time  to  approach 
90%  response.  So,  it  is  necessary  to  define  a  time  constant  tc  that 


Table  1 

The  values  of  the  parameters  used  in  this  work 


Physical  quantity 

Notation 

Unit 

Value 

Reference 

Viscosity  of  gas 

Pg 

kgm-1  s 

2.03  x  10"5 

[17] 

02  diffusivity  with  substance  of  Aira 

Do2 

2  -l 
m  s 

2.71  x  10“5 

Incropera  and  DeWitt,  p.927  [18] 

H2O  diffusivity  with  substance  of  Aira 

Dh2 

2  -l 
m  s 

3.35  x  10“5 

Incropera  and  DeWitt,  p.927  [18] 

Absolute  Permeability  in  CL 

kch 

2 

m 

1.2  x  10"15 

[17] 

Porosity  of  CL 

£CL 

0.4 

[17] 

Reaction  rate  constant 

K 

kmol  m-3  s 

1.5  x  10"5 

[11] 

Tortuosity  of  GDL/CL 

T 

1.5 

[11] 

Net  water  transport  coefficient  of  the  membrane 

a' 

0.5 

Set 

Universal  gas  constant 

R 

J  kmol- 1  K- 1 

8314 

Constant 

Faraday  constant 

F 

ckmol-1 

96487000 

Constant 

Operation  temperature 

T 

K 

353 

Set 

Transfer  coefficient  for  the  reaction 

a 

1 

Set 

Number  of  electrons  in  relation  participating 

n 

2 

Constant 

Thickness  of  GDL 

sd 

m 

3.0  x  10“4 

Set 

Thickness  of  CL 

Sc 

m 

1.5  x  10“5 

Set 

Thickness  of  membrane 

Sm 

m 

2.0  x  10“4 

Set 

Width  of  channel 

k 

m 

1.0  x  10“3 

Set 

Width  of  shoulder 

h 

m 

1.0  x  10“3 

Set 

a  The  values  are  derived  at  the  pressure  1  atm  and  temperature  353  K. 
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Fig.  2.  Effects  of  GDL  porosity  e\  on  time  evolution  of  average  current  density 
when  the  PEM  fuel  cell  is  subjected  to  a  start-up  operation  at  A P  =  7  mbar, 
7]  =  0.2  V. 

reflects  how  long  a  start-up  operation  or  a  state-to- state  opera¬ 
tion  will  almost  reach.  In  the  present  study,  the  time  constant 
is  defined  as  the  time  interval  needed  to  reach  90%  response  of 
average  current  density.  For  the  start-up  operation, 

i(tc )  =  0.9  i(t  =oc)  (19a) 


For  the  state-to-state  condition, 

i(tc)  =  O.l/i  +0.9/2  (19b) 

Where  i\,  ii  are  respectively,  the  average  current  density  of  the 
former  steady  state  and  the  latter  steady  state.  Start-up  operation 
can  be  considered  as  a  special  case  of  state-to-state  operation 
when  i\  is  zero. 

4.1.  Start-up  operation 

In  start-up  operation,  the  initial  current  density  is  zero  as 
initial  oxygen  mass  fraction  in  catalyst  layer  is  zero.  After  the 
system  is  started  up,  oxygen  mass  fraction  and  current  density 
gradually  rise,  and  soon  approach  a  certain  condition.  The  sys¬ 
tem  is  regarded  to  reach  a  steady  state  when  it  satisfies  a  certain 
criteria. 

Fig.  2  denotes  the  effects  of  GDL  porosity  £\  on  the  time 
evolution  of  average  current  density  when  the  PEM  fuel  cell  is 
subjected  to  a  start-up  operation  under  the  condition  of  pressure 
drop  7  mbar  and  surface  overpotential  0.2  V.  It  is  clear  that  as 
the  £\  increases,  the  average  current  density  in  steady  state  rises. 
This  means  that  better  mass  transport  in  the  PEM  fuel  cell  is 
noticed  for  a  larger  GDL  porosity.  In  addition,  the  response  time 
of  start-up  operation  is  fast  for  a  larger  GDL  porosity.  So,  a 


Fig.  3.  Evolution  of  oxygen  distribution  along  interface  between  CL  and  GDL  at  A P  =  1  mbar  and  r\  -  0.2  V  (a)  s\  =  0.3;  (b)  e\  -  0.4;  (c)  s\  =  0.6;  (d)  s\  =  0.8. 
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Fig.  4.  Evolutions  of  average  current  density  with  various  surface  overpotentials 
r\  at  A P  =  7  mbar  and  e\  =  0.6. 

higher  porosity  of  GDL  leads  to  better  cell  performance  at  the 
only  consideration  of  mass  transport.  The  value  of  time  constant 
is  quite  less,  about  0.1  s,  and  is  generally  less  than  that  of  the 
PEM  fuel  cell  with  straight  flow  field  [11]. 

As  the  catalyst  layer  is  quite  thin,  the  oxygen  mass  distri¬ 
bution  in  this  layer  can  be  regarded  one-dimensional,  and  the 
distribution  along  the  interface  between  the  catalyst  layer  and 


the  GDL  can  represents  the  distribution  in  whole  CL  region. 
Lig.  3  shows  the  evolutions  of  oxygen  distributions  along  the 
interface  between  CL  and  GDL  at  AP  =  7mbar  and  77  =  0.2  V 
under  different  values  of  s  1,  0.3,  0.4,  0.6  and  0.8.  The  time 
needed  for  almost  reaching  steady  state  from  start-up  condition 
can  be  figured  out  to  be  about  0. 1  s.  When  the  porosity  is  larger, 
uniformity  of  current  distribution  is  better. 

The  evolutions  of  average  current  density  with  various  sur¬ 
face  overpotentials  r\  under  a  certain  pressure  drop  7  mbar  and 
porosity  of  GDL  0.6  are  presented  in  Lig.  4.  A  careful  inspection 
of  Lig.  4  discloses  that  as  the  r\  increases,  the  average  current 
density  in  steady  state  rises,  and  the  time  constant  of  start¬ 
up  operation  decreases.  Additionally,  the  time  constant  varies 
largely,  ranging  from  0.01  to  0.1  s. 

The  evolutions  of  oxygen  distribution  along  the  interface 
between  CL  and  GDL  under  different  values  of  77,  0.175,  0.200, 
0.225  and  0.250  V  at  constant  pressure  drop  of  7  mbar  and  poros¬ 
ity  of  GDL  0.6  are  shown  in  Lig.  5.  It  is  clearly  found  in  Lig.  5  that 
the  times  needed  for  almost  reaching  steady  state  from  start-up 
for  different  surface  overpotentials  are  very  fast  and  range  from 
0.01  s  to  0.1  s. 

Ligs.  6  and  7  compare  the  effects  of  pressure  drop  on  the 
dynamical  behaviors  of  the  PEM  fuel  cells.  These  two  figures 
depict  that  the  steady- state  average  current  density  increases 


j  i  i - 1 - T"--  1  1  l  _ .  ■  ■ 

0.00  0.25  0.50  0.75  1.00 


Fig.  5.  Evolution  of  oxygen  distribution  along  interface  between  CL  and  GDL  at  AP  =  7mbar  and  £1  =0.3  (a)  77  =  0. 175  V;  (b)  ij  =  0.200  V;  (c)  rj  =  0.225  V;  (d) 
rj  =  0.250  V. 
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Fig.  6.  Effects  of  pressure  drop  on  the  dynamical  behaviors  of  the  current  density 
at  7]  =  0.20  V  and  e\  =  0.6. 

with  the  pressure  drop  and  the  response  time  decreases  with  pres¬ 
sure  drop.  The  difference  under  various  pressure  drops,  however, 
is  smaller  than  that  under  various  surface  overpotentials. 

Fig.  8  plots  and  compares  the  effects  of  the  surface  over¬ 
potential  o  and  the  pressure  drop  A P  on  the  time  constant  tc. 
Generally,  the  tc  decreases  when  the  ij  increases  or  when  A P 
increases,  and  the  value  of  tc  is  below  0.1  s.  When  the  value  of  rj 
is  relatively  small,  the  tc  changes  dramatically  with  A P,  and  the 
value  of  tc  at  A P  =  0  mbar  is  four  times  of  that  at  A P  =  10  mbar. 


Fig.  8.  Time  constant  variations  with  various  pressure  drops  A P  and  surface 
overpotential  at  e\  =  0.6. 

When  the  ij  is  relatively  large,  for  instance  larger  than  0.24  V, 
the  tc  changes  slightly,  and  remains  a  value  below  0.02  s. 

Although  the  initial  oxygen  mass  fraction  is  zero,  the  distri¬ 
butions  of  oxygen  change  at  once  after  the  system  is  started 
up.  The  oxygen  mass  is  always  not  uniformly  distributed  at 
the  beginning.  It  is  always  true  that  the  oxygen  mass  fraction 
near  the  inlet  is  larger  than  that  near  the  outlet.  The  deviations 
between  the  oxygen  mass  fractions  at  different  local  spots  vary 
with  time,  increasing  firstly  and  then  decreasing.  When  the  ij 
is  relatively  small,  the  oxygen  mass  distributions  develop  to  be 


Fig.  7.  Evolution  of  oxygen  distribution  along  interface  between  CF  and  GDF  at  rj  =  0.2  V  and  e\  =0.6  (a)  AP  =  0mbar;  (b)  AP  =  3mbar;  (c)  AP  =  7mbar;  (d) 
A  P  =10  mbar. 
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Table  2 


Corresponding  current  density  and  surface  overpotential  for  four  steady  states 


State  no. 

State  1 

State  2 

State  3 

State  4 

Current  density  (A/m2) 

4500 

10000 

15227 

20430 

Surface  overpotential  (V) 

0.206888 

0.22513 

0.239513 

0.256604 

quite  uniform.  When  the  r\  is  relatively  large,  the  distributions 
are  still  non-uniform.  It  seems  useful  to  introduce  the  ratio  of 
ORR  speed  to  gas  transport  speed  to  account  for  the  discrepancy. 
The  larger  the  value  of  the  ratio  is,  the  more  non-uniform  the 
distribution  will  be.  In  the  corresponding  conditions  of  Fig.  5, 
the  pressure  drop  is  constant,  which  means  similar  gas  trans¬ 
port  process,  and  r\  dominates  ORR  speed  so  as  to  dominate  the 
ratio.  In  the  corresponding  conditions  of  Fig.  7,  surface  over¬ 
potential  is  constant,  meaning  similar  ORR  speed,  and  pressure 
drop  dominates  transport  speed  so  as  to  dominate  the  ratio. 

4.2.  State -to- state  operation 

In  the  state- to- state  operation,  the  initial  conditions  differ 
from  that  of  start-up  operation.  The  main  feature  is  the  oxy¬ 
gen  mass  distribution  is  not  void  but  at  a  steady  condition  of  the 
former  state.  In  practical  operation,  operators  change  operating 
condition  through  external  load  adjustment,  and  the  overpoten¬ 
tial  reflects  the  current  drawn.  So  the  current  density,  in  this 
study,  is  chosen  as  the  identification  variant  of  various  states. 
Stable  cases  analysis  is  firstly  conducted.  Fig.  9  shows  the  polar¬ 
ization  curves  for  various  pressure  drops  A P.  It  indicates  that  the 
limiting  current  density  increases  with  the  pressure  drop.  It  is 
quite  normal,  as  larger  pressure  drop  means  better  gas  trans¬ 
port  capability  and  less  transportation  overpotential.  Fig.  10 
expresses  the  relation  between  surface  overpotential  and  cur¬ 
rent  density  for  various  A P.  According  to  the  plot  of  Fig.  10, 
four  states  at  A P  =  7  mbar  are  chosen  to  analyze  transient  char¬ 
acteristics  of  state-to- state  operation.  The  detailed  information 
about  the  four  states  is  listed  in  Table  2. 

Fig.  1 1  shows  the  dynamical  performance  from  one  state  to 
another  state.  When  the  cell  performance  alters  from  a  lower 
current  density  state  to  a  higher  current  density  state,  as  the 
upper  graph  depicts,  the  current  density  jumps  up  to  a  quite  high 
level  in  the  beginning,  and  then  decreases  gradually  till  to  the 


Table  3 


Response  time  constants  from  state-to-state  operation  (Unit:  s) 


start-up 

State  1 

State  2 

State  3 

State  4 

start-up 

0.024803 

0.015634 

0.00889 

0.004333 

State  1 

0.016167 

0.009272 

0.004230 

State  2 

0.025046 

0.009109 

0.004084 

State  3 

0.024463 

0.015564 

0.003879 

State  4 

0.023891 

0.015107 

0.008432 

steady  state.  On  the  contrary,  when  the  cell  performance  alters 
from  a  higher  current  state  to  a  lower  state,  as  the  lower  graph 
depicts,  the  current  density  jumps  down  to  a  quite  low  level  in 
the  beginning,  and  then  increase  gradually.  For  both  operations, 
the  value  of  current  density  may  change  up  to  5  times. 


Fig.  10.  Relation  between  current  density  and  surface  overpotential  for  various 
pressure  drops  A P. 


Time  [s] 


Fig.  11.  Dynamic  cell  performance  from  one  state  to  another  state. 
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The  time  constants  from  one  state  to  another  state  are  listed 
in  Table  3.  The  values  of  time  constants  also  differ  obviously. 
The  large  value  is  one  order  higher  than  the  low  one.  As  shown 
in  Table  3,  the  response  time  generally  depends  on  the  new  state. 
For  instance,  the  response  time  is  around  0.024  s  when  the  new 
state  is  state  1 .  It  is  necessary  to  notice  that  the  time  constant 
from  state  A  to  state  B  is  different  with  that  from  state  B  to  state 
A.  For  example,  the  time  constant  tc  from  state  1  to  state  4  is 
0.004230  s,  and  tc  from  state  4  to  state  1  is  0.023891  s. 

5.  Conclusions 

In  this  work,  the  transient  transport  process  of  PEM  fuel  cell 
with  interdigitated  flow  field  is  numerically  investigated.  Two 
operations,  start-up  and  state- to- state,  are  analyzed.  Parametric 
study  is  performed.  Some  conclusions  are  resulted. 

1 .  The  transient  response  time  is  generally  quite  small,  below 
0. 1  s,  and  decreases  with  the  surface  overpotential,  the  poros¬ 
ity  of  GDL  and  the  pressure  drop.  But  when  the  surface  over¬ 
potential  is  considerably  large,  the  response  time  changes 
slightly  with  the  pressure  drop. 

2.  The  ORR  speed  and  mass  transport  speed  determine  the 
response  time.  Response  time  of  the  PEM  fuel  cell  with 
interdigitated  flow  field  is  quite  less  than  that  of  the  PEM  fuel 
cell  with  straight  flow  field,  as  the  mass  transport  speed  of  the 
former  cell  is  quite  larger  than  that  of  the  latter  one.  The  mass 
transport  of  the  former  cell  is  dominated  by  both  convection 
and  diffusion,  and  that  of  the  latter  is  only  by  diffusion. 

3.  The  ratio  of  ORR  speed  to  mass  transport  speed  determines 
the  steady  oxygen  mass  distribution  in  catalyst  layer.  The 
larger  the  value  of  the  ratio  is,  the  more  non-uniform  the  dis¬ 
tribution  will  be.  A  relatively  low  overpotential,  a  high  poros¬ 
ity  of  GDL  or  a  high  pressure  drop  leads  to  a  better  uniformity 
in  current  density  distribution  along  the  width  of  the  cell. 

4.  For  state- to- state  operation,  the  value  of  current  density 
changes  largely,  especially  in  the  beginning.  The  response 
time  mainly  depends  on  the  new  state  condition. 
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